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Abstract 

We study how the encasement of a growing elastic bulk within a possibly differ¬ 
ently growing elastic coat may induce mechanical instabilities in the equilibrium 
shape of the combined body. The inhomogeneities induced in an incompressible 
bulk during growth are also discussed. These effects are illustrated through a 
simple example in which a growing elastic cylinder may undergo a shape tran¬ 
sition towards a bent configuration. 
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1. Introduction 

Shape is one of the most prominent features of cells, plants, organs, and 
of living beings in general. Specific shapes may play a crucial role in enabling 
the achievement of specialized biological goals, which might allow for an organ¬ 
ism’s survival and thriving in its environment. Morphogenesis has been linked 
to mechanical origins in a natural way because mechanics offers innumerable 
examples of systems whose equilibrium configurations depend on external pa¬ 
rameters, with bifurcations among equilibria which may imply very different 
shapes and structures for a body. Perhaps the most studied and exploited ex¬ 
ample of such morphogenetic bifurcation phenomena is given by Euler’s elastica. 
In more recent times we mention the studies by Biot [5] who understood how 
three-dimensional rubbers under compression may develop surface instabilities, 
as have also later been found in soft, strain-hardening materials [71 Hi]. 
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A mechanical approach in the creation of form for bodies undergoing growth 
was envisaged already at the beginning of last century by Thompson [SS]. In¬ 
teresting studies along these lines were later performed, such as for instance [9], 
where shell theory was used to estimate the surface stresses that can be induced 
in a growing surface, showing how they may related to splitting and cracking in 
vegetables. These methods give very fertile ground for the mechanical investiga¬ 
tion of bio-morphogenesis, greatly expanding the scope of continuum mechanics 
well into the boundaries of biophysics and biology in general. A growing body 
indeed evolves in a natural way under the effect of a number of external pa¬ 
rameters which vary with time inducing morphogenetic effects with lesser need 
for genetic encoding of information. This offers attractive avenues for the ex¬ 
planation of form ISO], within a background of natural hypotheses of genetic 
parsimony. 

Understanding in terms of mechanical instabilities has since been proposed 
for numerous observed shapes of growing bodies, using a wide variety of ap¬ 
proaches, and sometimes exhibiting remarkable agreement with experimental 
observations. To name but a few examples, buckling under external stress has 
been related to morphogenesis for instance in growing spheroidal shells m , or 
in the growth of other constrained systems dalle]. More complex shapes, ob¬ 
served in plant organs such as long leaves [53] ? or in blooming flowers mm, 
have also been modelled by adapting the theory of elastic shells to a growing sur¬ 
face. See da for a review of the mechanics of buckling-related morphogenesis. 
Further recent studies have further analyzed the growth of vegetable matter, as 
in growing pumpkins [21] or in the ripening of kiwi fruit [18] , see also the review 
[28] . In the biomedical sciences, the investigation of mechanical instabilities has 
helped investigate the role of external stress in tumor growth mi, and me¬ 
chanics has since provided fruitful soil for cancer modelling [29] . Recent work 
has also related mechanical instabilities to morphogenetic development during 
the growth of animal organs [331136] , and in brain formation [i[20l[32]. 

In most cases, the basic trigger for morphogenetic instabilities derives from 
the presence of differential growth and distinct mechanical properties in different 
neighboring parts of a growing body. Thus, adjacent layers or domains sharing 
open two-dimensional mm or one-dimensional [33] adhesion boundaries have 
been examined. Confinement of a growing body within a closed regular two- 
dimensional boundary not undergoing growth has also been considered [241131] . 
Less explored appear to be the morphogenetic possibilities for encased growing 
bodies, wherein a growing bulk is surrounded/protected by a closed regular 
growing surface layer (as with the skin or rind of fruits), whose mechanical and 
growth properties will in general differ from those in the bulk. 

The aim of the present study is thus to highlight how encasement may gener¬ 
ate interesting morphogenetic instabilities in growing bodies. We treat growth 
by following the approach introduced in mm, and extensively adopted in 
theoretical studies of growing systems, as in gidn] among many examples. A 
peculiarity in our approach is that two multiplicative decompositions for the 
deformation gradient must be introduced to correctly account for the possibly 
independent growth of the bulk and the coat. In the bulk, growth is described 
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by the factor Fg in the decomposition F = F^Fg of the three-dimensional de¬ 
formation gradient F. Simultaneously, a similar decomposition A = AgAg 
is introduced for the surface deformation gradient A. The strain energy then 
depends only on the elastic factors Fg and Ag. By following [niiss] in our 
treatment of elastic continua surrounded by elastic surfaces, we consider an in¬ 
compressible neo-Hookean material for the bulk, and Foppl-Van Karman elastic 
coat, undergoing finite strains. 

The growth components Fg, Ag are at least partially determined by external 
processes, typically of biochemical origin. Furthermore, material remodeling 
may possibly influence such processes, and at least partially relieve the residual 
stresses which may have been originated by incompatible growth. The choice 
of including or not such an evolution for the growth components determines 
whether we are considering morphogenesis at stress-free conditions vs. in the 
presence of residual stresses. This might for instance depend on the relative time 
scales of growth and relaxation (see [5] for an example of how a soft material 
may exhibit purely elastic or apparently plastic behavior, depending on the 
deformation and relaxation times). We keep our study as simple as possible 
here, and do not consider any remodeling, so that in general the deformed 
shapes we characterize possess residual stresses (as, for instance, the systems 
considered by [H]). 

Our analysis is based on the combination of two effects. The first is the role 
that non-homogeneous growth may play in giving incompressible bodies the 
possibility of undergoing apparently non-isochoric deformations. The second 
factor is showing how encasement, i.e. the coexistence of a growing bulk inside 
a possibly differently growing skin, may affect the equilibrium configurations 
of the body even in the absence of instabilities of the bulk or the boundary 
separately. 

Here we illustrate these effects by means of a simplest example, that is, a 
bending instability for an isotropically growing straight cylinder encased within 
a cylindrical boundary which is in turn growing at possibly different rates in 
the longitudinal and transverse directions. Such growth anisotropy may occur, 
for instance, when parallel fibers are present in the skin, which may also possess 
anisotropic elastic moduli. We show how encasement induces, in a range of 
growth and elastic parameters, instabilities which lead to new equilibria with a 
bent shape for the cylinder. Depending on the material properties, other insta¬ 
bilities, which we do not explore here, may occur in encased growing cylinders, 
such as bulging, twisting, or others (see for instance the deformations consid¬ 
ered in [TS]). Our analysis shows that for a growing body the constraint of 
encasement is different from simple confinement as in EllSI], as the growth of 
the external coat along with the growth of the bulk creates instabilities which 
might not be present when the coat is (deformable but) incapable of growth. 

This paper is organized as follows. In Sectionj^we show how inhomogeneous 
growth may provide new degrees of freedom for the distortion of incompressible 
growing bodies. In Section we analyze the effects of encasement on growth, 
showing how mechanical instabilities may arise when a growing body is sur¬ 
rounded by a differently growing coat. In Section we evidence a bending 
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instability generated by encasement. An Appendix contains some of the com¬ 
putations leading to the results presented in the main text. 

2. Growth in incompressible materials 

In this section we analyze how inhomogeneous growth may enlarge the class 
of deformations available to incompressible materials. Let x be the map trans¬ 
forming the reference configuration into the present placement, and = Vy. 
We model growth by assuming [Hill] that F may be decomposed as the prod¬ 
uct of an elastic and a growth component: 

F = FeFg. (1) 

The growth tensor Fg describes how growth would locally shape body volume 
elements, were they allowed to grow stress-free, that is, in the absence of the 
surrounding elements. The balance equations driving Fg may not entirely be 
included in the bio-mechanical model. When this is the case, some or all of its 
entries may be explicitly specified, and therefore treated as external parameters. 
We remark that Fg need not be the gradient of any ‘growth deformation’, as 
growth is a local phenomenon which may induce a lack of compatibility when¬ 
ever curling ^ 0. Moreover, since mass is not conserved during growth, no 
isochoricity condition must be a priori enforced on deti^g. 

The elastic distortion Ff. is defined by the decomposition 0: Fe = FF-\ 
The assumption that Fg identifies the current stress-free configuration implies 
that the strain-energy density must be a function of the elastic strain 

Ce = Fj Fe = Fg^F^FFg^ = F~^CF-\ (2) 

where C = F^F is the standard Lagrangian strain tensor. As Ce is itself a 
function of G, any strain-energy density depending on Ce automatically com¬ 
plies with frame-invariance requirements. 

The incompressibility constraint establishes an isochoricity condition on the 
elastic distortion Fe, so that 


detF = detFg. (3) 

If the point-wise mass supply is assigned, the map det Fg must be treated as 
an external parameter, and incompressibility limits as usual the class of defor¬ 
mations available to the growing body. If, on the contrary, the system deforms 
sufficiently slowly, apparently non-isochoric deformations become possible, in 
which growth accommodate local volume variations by concentrating mass pro¬ 
duction at sites undergoing greater expansion detFg. When this is the case, 
the incompressibility constraint applies only globally, and the determinant of 
F may neither be 1, nor even uniform. Much more general deformations thus 
become available to a slowly growing body, the map det F giving the explicit 
details on the local mass production. 
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Figure 1: (Left) Reference straight capped encased cylinder. (Right) Grown, stretched and 
bent configuration for the same cylinder. 


3. Encasement 

We now study how the independent growth of an elastic bulk and of a closed, 
regular elastic surface surrounding it, may induce morphological instabilities. 
We consider the simplest geometry and constitutive assumptions on the material 
to better evidence the origin of such instabilities. The bulk of the body is given 
by a capped, isotropically growing cylinder made of simple incompressible neo- 
Hookean material, coated with a possibly anisotropic Fdppl-von Karman elastic 
surface undergoing finite strains. Two half spherical caps are considered, with 
centers on the cylinder axis and suitable radius, to enforce encasement, while 
ensuring C^-regularity of the external boundary, see Fig. left. We denote 
rg, /ig respectively the cylinder radius and height in the reference configuration, 
which in cylindrical coordinates is given by 

Bo = {P € £ ■. P = O + pcos 0 63 ; + psinO Ey + Z e^, 

Z G (—ro, ho + rg), p G [O, Pro,ho (Z )), 0 € [0, Stt)}, 

where {ex, is a fixed orthogonal basis and, for any positive r 

Rr,h ■ {—r, h + r) ^ M+ is defined as 

rVr 2 -z 2 if 2 ;g(-r, 0 ] 

Rr,h{z) = < r if z e ( 0 , /i) 

[\/r2 — (z — /i)2 \i z G\h,h + r). 

The reference body is therefore identified by the reference radius rg and height 
ho of its cylindrical part. 

As mentioned in the Introduction, we are interested in determining under 
which conditions the straight shapes of the capped cylinder become unstable, in 
particular in favor of bent stable shapes. We perform our analysis by introducing 
a parameterized family of shapes, and checking the stability of straight cylinder 


(4) 

and /i, 

(5) 
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configurations against bending. The critical values we identify in this way thus 
provide a sufficient criterion for the loss of stability of the straight configuration. 
Other instabilities, not presently analyzed, could in principle occur at different 
critical thresholds, depending on the values of the constitutive parameters. 

In detail, the family of bent shapes that we consider is obtained by assuming 
that the axis of the cylinder becomes the arc of a circle with curvature k, as in 
the right panel of Fig. (see the Appendix for further details): 

Bk. = {P € £ '. P = Qk + pcosOe^ + (psind — (cos kz — sinreze^), 
z e (-r, h + r), pS [0, i?r,h(z)), 9 G [0, 27r)}, (6) 

where Qk, = O + kT^ By is the center of the circle on which the axis lies. The 
domain Bk, approaches the upright shape Bq for vanishing k. 

The encased cylinder undergoes growth, with details on the bulk and sur¬ 
face growth distortions Fg, Ag provided below, while it deforms to achieve the 
current configuration B,^, determined by minimizing the strain energy below. 

• The bulk material undergoes a possibly inhomogeneous dilation, charac¬ 
terized by the growth tensor (see (§) 

^g|buik = (detF)^/"/, (7) 

where I is the identity. Since the bulk material is assumed to be in¬ 
compressible, growth fixes the total volume of the current configuration . 
We therefore introduced a positive scalar a which represents the relative 
volume increase in the grown body, and require 

Kurr = (8) 

• The coat of the cylinder undergoes an independent growth process, de¬ 
scribed by anisotropic elongation factors /?, 7 , respectively in the longi¬ 
tudinal and radial directions in the cylindrical part of the coat, plus an 
isotropic growth factor 7 for the semi-spherical caps. We do not require 
the coat material to be inextensible, although stretches with respect to 
the locally grown configuration will be energetically penalized. 

In this framework, growth is therefore encoded in the geometrical dilation 
factors a, /3, 7 , fixed by external processes, which we treat as parameters in the 
analyzed solutions of the problem. We remark that encasement induces a lack 
of compatibility in the field Fg whenever (3 differs from 7 , as no dilation of the 
growing bulk may succeed in accommodating anisotropic surface growth. This 
occurs often in growth processes, and is at the basis of the existence of residual 
stresses [32] . 

3.1. Strain energy 

To determine the optimal deformed configuration, we follow [35] and assume 
that the total strain energy of the encased cylinder can be decomposed into a 
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bulk and a surface contribution 


J^[X]= f a{C)JgdV+ [ ,{A,K)ZsdA, 
JBo JdBo 


(9) 


in the absence of any loads. The strain-energy density depends in fact only 
on the elastic distortion but, by using ([^ for the bulk term and a similar de¬ 
composition for the surface deformation gradient (see (A.20) in the Appendix), 
it may be written as density in the reference configuration, depending on the 
Lagrangian strain C, the surface deformation gradient A, and the curvature 
tensor K. The Jacobians Jg,5g, respectively of the bulk and surface growth 
transformations are included in ([^ because the strain-energy density is defined 
with the respect to grown volume/area elements. 

We assume a simple, isotropic neo-Hookean behavior for the bulk strain- 
energy density: 


a(C) =^(trCe-3) =At(tr - s) 

= ^^tr ((det — 3^ = /r(trC — 3) 


( 10 ) 


with ^ the shear modulus, and C = (det the isochoric strain. 

To include in our study materials with possibly anisotropic (e.g., fiber- 
structured) coats [34], we choose the following anisotropic expression for the 
surface strain-energy density 


c(A, K) = ,S|| {E^T^ ■ T^f + Sa_{E,T^ ■ + BZeH 


( 11 ) 


where ,T^ are respectively the longitudinal and azimuthal tangent unit vec¬ 
tors on the reference surface (see Appendix A.2|), and Sj_ represent the lon¬ 
gitudinal/transverse 2D stiffness moduli, and B the bending stiffness. The first 
and second terms in the surface strain-energy 0 depend on the elastic com¬ 
ponent Ae of the surface deformation gradient A = A^Ag, as Eg = AjA^ — Is, 
where Ig denotes the identity operator on the tangent planes to Bq. The term 
penalizing the mean curvature H carries a further factor Ze = detAo, as it 
penalizes curvature in the current configuration. Eq. 0 defines for the coat a 
finite-elasticity counterpart of the strain energy of a classical Foppl-von Karman 
plate |39|. 


4. Bending instability 

We now examine under which conditions the straight configuration may be¬ 
come unstable, possibly leading the cylinder to undergo a morphological trans¬ 
formation towards a bent shape (some of the computations are given in the 
Appendix). We remark that our model treats separately the bulk and surface 
strains, so that fairly straightforward changes of the computations presented 
below could also account for boundary slip and/or detachment during encased- 
growth phenomena. However, these are absent in the ensuing first analysis of 
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encasement instabilities, as the family of deformations which we consider now 
is continuous at the coat-bulk boundary. 

We decompose the capped cylinder into the cylindrical part, plus the two 
spherical caps, and consider the family of deformations x explicitly defined in 
the Appendix (see (A.4|), parameterized by the height and radius h,r of the 
current configuration, and the curvature k of the cylinder axis. The deformation 
gradient of the bulk transformation in the cylindrical part of the material (see 
(A.4)) may be written in the orthogonal basis used in Q as 


FiX,Y,Z) = -(e, 
ro 

with 


f{Z)®ey) 


1 - 


nrY 

ro 




nhZ . kHZ 
j{Z) = cos —-— By — sm —-— fiz 


. KhZ 

f±[Z) = sm-—ey 


Khz 


cos 


( 12 ) 


(13) 


The Jacobian 


J = det F = 


hr^{rQ — ktY) 
horl 


(14) 


provides the volumetric dilation factor in the bulk, and sets a maximum admis¬ 
sible value for the curvature k for the deformation to be regular. More precisely, 
since |T| < tq, the deformation is regular only when Kr < 1, so that 

0 < K < r~^. (15) 

The cylindrical contribution to the volumetric strain-energy density gives 
(see (10 1 and (A.9)) 

I ^hr^{ro — ktY) 


Icyl 


^ 0^0 


/ hpr \ i / h{ro - KrY) \ I _ 
\h{ro — KrY)) V hnr ) 


hor 


(16) 
= I, and 


The caps only undergo a dilation by a factor r/ro, so that C\ 
a\ = 0 . 

I caps 

Identifying the elastic component of the surface deformation gradient 
and the mean curvature H requires lengthier computations, reported in the 
Appendix. In the deformed conhguration, the cylindrical section of the coat 
stores a strain energy given by (see (A.20 1 and (A. 21 |): 

2 


“’Icyl 4 


^ii / h‘^{l — Krsin0)^ 




- 1 




7^2 


- 1 


Bh 


(17) 


4horor(l — wrsin©) 

The coat caps undergo an isotropic extension of factor r/ro, remaining half¬ 
spheres of mean curvature H = —r~^ (see (A.23l). Therefore 


(^11 +S±) f r 




-1) +4 


2 ■ 


(18) 

























Let us now introduce the dimensionless parameters u = r/(aro), v = h/{aho), 
which reflect whether the current shape respects the dimensions obtained by a 
simple rescaling of the reference shape with a factor dictated by Q . Let further 
K = utqk measure the bending state of the current configuration, and (5 = j5/a, 
7 = 7/0 be the rescaled surface growth factors. We also introduce the reference 
aspect ratio Aq = ro/hg. Notice that (u, v, k) parameterize the current con¬ 
figuration, while the parameters (/?, 7 ) characterize the normalized anisotropic 
growth. The total strain energy takes the form 


J^= [ a{F,)J^dV+ f ^{A,K)Zg,dA 
JBo JdBo 


I Bo 
,3i, „2„,2 


= fia'^hor^u V / pdp 


de 


2u 3 (1 — Ku sin 0) 3 u3 (1 — KMsin 0) 3 


-I- a^roho/3^ / dQ 


S'!! /u^(l — Kusin 0 ) 


V 3 

PL 2 


U3 

2 ^ 2 -, 


-II 


Bvho 

Auro 


r2'w 


de 


1 — ku sin 0 


-|-7rQ;"'rgy(S'||-|-S'_L) +ATrB. 


To compute the above integrals we define, for any |a| < 1 , 

fi{a,b) = ^ J (1 - asin0)^d0 = 2F1 (5(1 - b ), -56; 1; a^) 

with 2T1 the Gauss hypergeometric function such that 

/i(a,0) =/i(a, 1) = 1, and 

/i(a, 6) = 1 -I- — 1 ) -I- 0 (a"^) as a —>■ 0 . 

For future use, we also introduce, for any |a| < 1 , 

f2(a,b) = 2 f pfi{ap,b)dp = 2F1 (^(1 - b) , -^b; 2; a^) . 
^0 


We thus obtain 


T — 21^ ixhi^r^o? li^ V 


ri r ^ 


/o 


2^fiikup, i) -f ^/i(kup, |) - 3 

V3 U3 


dp 


+ 27ra^ro/io/37 
TTBvhQfi{ku, — 1 ) 


/.,4 

4 


+ 


2u^/i( kM, 2) \ S± f 


u 


(19) 


-3 


( 20 ) 


( 21 ) 


( 22 ) 


- 1 


T 


2urQ 


-HdTrS. 


7 


(23) 


The hypergeometric integrals in (23) may again be computed, and provide the 
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strain energy per unit current volume ('I' = F/(Tra^rQho)) 

'I' = ^(2uSn3/2(KM, |) + u^vi f2{iiu, |) — Sm^u) 

, 2/37 


aro 


S\\ (v'^fi{k.u,4:) _ 2v‘^fi{ku,2) ^ 5'_l ^ 


/34 


2 1 


Bvfi{Ku,-l) Tecs'll + S'_l)Ao 


2 argM 


aro 


- 1 


4BAn 


3^-3 ■ 


a^’r. 


(24) 


We remind that u and v are not independent variables. Indeed, the volume 
of the current configuration is fixed by the growth, so that (see equation 


u^{v + |Aou) = 1 + |Ao, 
where Aq represents the reference aspect ratio. 


(25) 


Eq. (241 evidences that encasement, i.e. the combined effect of bulk and coat 


growth, may not be simply renormalized as the effective growth of one of them 
(either bulk or coat) with the other considered as a simple elastic medium, so 
that encasement does not in general reduce to the constraint of confinement. 
For instance, in the present case, we describe the surface growth in terms of the 
renormalized parameters /3 and 7 , but the bulk growth parameter a enters also 
the ratios between surface and bulk elastic constants. Therefore, when a (i.e. 
the bulk) grows, even if the relative growth parameters /3 and 7 remain constant, 
the surface effects become less important, affecting the stability thresholds and 
possibly their existence. 

The straight shape is always a stationary configuration for the energy func¬ 
tional, because, setting d'o = 


4* = 4'o -b 


fiu3v3 [7v‘^ — u'^) S'||U^?;^(3n^ — , 3^)7 Buv 


18 


2arQl3^ 


+ 


Aarg 


k? + o(k^) (26) 


as if —> 0 , with 


'2\2 


dto = + 2u3v 3 — 3u^v) + 


S'_L/3(u2 - 7^)2 5'||^(^2_^2) 


2 aro 7 ^ 


2aro/3^ 


By 

2arQU ho \ 2aro^'^ 


2B 

aVg 


(27) 


To check the stability of the straight shape we assume that the unknown 
curvature k is small enough that u may be determined by minimizing the zero- 


curvature strain energy dtg, as defined in (271. We then replace the values 
obtained in the square-brackets term in ( 261 , checking its sign to test stability. 


This leads to a number of possibilities for the bending instability of the growing 
encased cylinder, depending on the imposed material parameters. 

Among these possibilities, we only examine here a particular case. We 
look for the equilibrium configurations of a cylinder with initial aspect ratio 
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unstable 



Figure 2: Phase diagram evidencing the onset of an instability for the straight con¬ 
figuration when the longitudinal growth, represented by /3, is large enough. Displayed 
values of the longitudinal-to-tangential stiffness ratio are, top to bottom, 5'||/S'j_ = 

0.1,0.15,0.2,0.25,0.3,0.35, with fixed 5 _l. 


Aq = To/ho = 1 , with a soft (though incompressible) bulk, almost negligible 
bending stiffness and an anisotropic coat, reinforced on the tangential direction. 
Precisely we set B/{2arQfi) = 10“^, S±/{2arQfi) = 10^, and vary as indi¬ 
cated in Fig.l2]For each value of and 7 we determine numerically the optimal 
value of u in p7|), and check the straight stability sign in (26). Fig. [^evidences 
how, for each 7 , a critical value of /3 exists, above which the straight configura¬ 
tion becomes unstable. Due to the meaning of the growth constants 7 and /3, 
this in turn indicates the critical value of the aspect ratio of the grown cylinder 
at which there is the onset of the bending instability. We notice in Fig. that 
such critical value increases as approaches S±, indicating that surface elastic 
anisotropy helps triggering the instability of the straight configuration. 

The features of the bifurcation connected with such instability may be ana¬ 


lyzed by using the energy functional (24), obtaining for each value of 7 and (3 


the optimal value of the axis curvature k. We use the same material parameters 
as in Fig. along with S'||/S'_l = 10“^, and 7 = 0.8. We see from Fig. [^that 
the bending transition is of the second order, with the curvature continuously 
departing from zero (the underlying bifurcation for the growing encased cylinder 
corresponds to an axially-symmetric supercritical pitchfork). 
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Figure 3: Optimal value of the curvature k as a function of /3 for the material parameters 
provided in the text. The critical value of the longitudinal growth is /3 = 5.617. The dashed 

line is a best fit for the critical behavior, showing that itopt ~ sj0 ~ as /3 approaches the 
critical value leading to the bending bifurcation. The same parameters as in Fig. are used. 
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Appendix A. Geometry of bent cylindrical deformations 

Appendix A.l. Bulk deformation 

Let us consider the reference and current illustrated in Fig. 

Bq = {P £ £ : P = O pcosOe^ psinOe^ -|- ze^, 

z G {-ro,ho -f To), p G [0,Rro,ho{z)), © G [0, 27r)}, (A.l) 

Bk. = {P & £ ■ P = Qk + pcosOe^, -I- (psinO — k~^) {cos kz By — siuKZBz), 
z G [-r, /i -t r], p G [0, ( 2 )] , o G [0, 27r)}, (A.2) 

where ro, hg, r, h, and k are positive parameters, O and = 0-1- k~^ By 
are points in the three-dimensional Euclidean space £ with orthogonal basis 
{bx, By, B^}, and Rr^h ■ {—r, hr) ^ K+ is defined as 


(^/r'^ — z'^ it z G {—r, 0] 

BrM{z)=<r itzG{0,h) (A.3) 

— (z —/i)2 ii z G [h,hr). 

In all our calculations we consider the smooth deformation x '■ ^0 ^ B^ 
such that, for every P = O -I- Y e^; -|- Y -|- Z G Pqj the coordinates of the 
transformed point x{P) = O x b^ y By z Bz are given by 



rX 


X = 

— 



To 



rY 

/ 1 rY 

y = 

-h 1 



ro 

V ro 


hZ 

/ , rY 

z = 

T-h 

(«"- 


ho 

ro 
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(A.4) 






We introduce the orthogonal unit vectors 


KhZ nhZ p / \ k/iZ nhZ 

—— By — sin —— Bz and jj_[Z) = sin —— By + cos —— 


f{Z) = cos 

ILQ ILQ ILQ no 

(A.5) 

which, along with constitute an orthogonal basis in £. The deformation 

s 

f_L{Z) (g) Bz 


gradient associated with (A.4) may thus be written as 

ktY \ 


F = —(e^ 
ro 


f{Z) ® By 


1 - 


ro 


-) 


The isochoric deformation gradient is defined by 


F = J 3 F with J = det F = 


hr^(rQ — KrY) 
hoTo 


and the related isochoric right Cauchy-Green strain C is 

hor ^ \ , (h(rQ — KrY)\^/^ 


C = F' F = 


h{rQ — KrY) / 


/ f 

{I- Bz® Bz) + (^- 


hor 


(A.6) 


(A.7) 


) Bz. 

(A.8) 


In particular, 


trC = 2 


hor 


h{ro — KrY). 


(a^9) 


We end this computation by determining the condition on the parameters 
/i, r which ensures that the current placement occupies a prescribed vol¬ 
ume I4urr = ct^Vo, where Vq = Cq(/io + |?’o) is the reference volume, and the 
positive scalar a provides the bulk volume increase. The current caps are two 
hemispheres of radius r, and so their volume is |7rr^. The volume of the bent 
cylindrical part in can be computed as 


r^o 


r2Tr 


i,j pip! 

Jo Jo hor^ 


(A.IO) 


and does therefore not depend on the curvature k. As a consequence, the volume 
grows by a factor of a if and only if 


r\h + |r) = a'^roiho -f |ro). 


(A.ll) 


Appendix A.2. Surface deformation 

We now compute the surface deformation gradient which maps tangent vec¬ 
tors on the reference coat dBo into tangent vectors on the boundary of the 
current deformed body dB^. We refer to |3S] for more details. We first consider 
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the deformation which maps the cylindrical part of dB^ into the corresponding 
portion of dB^.- We introduce the notations: 

6^.(0)= cos 0 e^; + sin 0 
ee(0) = — sin0 e^; + cos 0 

(A.12) 

h{Z) = cos (KhZ/ho)ey — sin {KhZ/ho)e^ 
h±{Z) = sin {KhZ/hQ)ey + cos (KhZ/ho)ez, 

and identify points on dBo by the coordinates (0,Z) G {[0,27r) x [0 ,/iq]}, with 
1^(0, Z) = O + ro 6^(0) + Z Bz- Then, y(0, Z) = x(Y(0, Z)) is given by 

y(0, Z) = Qi^ + rcosQe^ + (rsin0 — K~^)h{Z), (A.13) 

where again = O + k~^ Sy. The tangent unit vectors in dBo and dB^,, which 
in this simple example may be identified with their duals, are respectively given 

by 

Te = T® = ee, te = = — sin0 e^, + cos0/i(Z) 

Tz = = Bz, tz = = h±{Z), 

thus leading to the following grown and current normal unit vectors 

N = Br n = cos 0 e^; + sin 0 h{Z). 


(A.14) 

(A.15) 


The metrics associated with the parameterizations provide the local area dila¬ 
tion factor, which connects the current and the reference area element: da = 
5(0) dA, with 


3 ( 0 ) 


hr{l — nr sin 0) 
horo 


(A.16) 


The surface deformation gradient 


A = —tQ®T^ 
ro 


h{l — Kr sin 0) 

ho 


(A.17) 


may be decomposed into an elastic and a growth component: A = A^Ag. This 
latter represents an anisotropic elongation of factors /3, 7 , so that 

Ag= 7 T®(g)^®-h/ 3 ^^(g)^^ and (A.18) 

Ac = AA-^ = —te(g)T°-h-^- — -^t^(g)TA (A.19 

"fro pho 


The elastic strain tensor is then 


= KA 
_ 1 
“ 2 


Ae - Is) 


7^2 


-I T' 


ft,2(l — Kr sin 0)2 


-IT 


, (A.20) 


where Is = T® 0 T® is the identity tensor in the reference tangent 

plane. 
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We obtain the mean curvature H from the current tensor curvature 


iJ = tr (g) 


Ksin0 
1 — nr sin 0 



1 

2r(l — nr sin 0) 


(A.21) 


The deformation in the caps is simpler, as it transforms a reference half¬ 
sphere of radius tq into the current half-sphere of radius r. Since Agjcaps = lIs, 
we obtain da = Z dA, with 


pans — 2’ pans — and Eg 

I Co-pS C3-pS ^^^^0 


I caps 


7^2 


-I T 


(A.22) 


The curvature tensor is isotropic, with mean curvature 


H 

caps 


—r 


-1 


(A.23) 
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